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ABSTRACT 

We present the first near-IR scattered light detection of the transitional disk 
associated with the Herbig Ae star MWC 758 using data obtained as part 
of the Strategic Exploration of Exoplanets and Disks with Subaru, and l.l//m 
HST/NICMOS data. While sub- millimeter studies suggested there is a dust- 
depleted cavity with r=0'.'35, we find scattered light as close as 0.1" (20-28 AU) 
from the star, with no visible cavity at H, K', or K^. We find two small-scaled 
spiral structures which asymmetrically shadow the outer disk. We model one 
of the spirals using spiral density wave theory, and derive a disk aspect ratio 
of h~0.18, indicating a dynamically warm disk. If the spiral pattern is excited 
by a perturber, we estimate its mass to be 5^4 Mj, in the range where planet 
filtration models predict accretion continuing onto the star. Using a combination 
of non-redundant aperture masking data at L' and angular differential imaging 
with Locally Optimized Combination of Images at K' and Kg, we exclude stel- 
lar or massive brown dwarf companions within 300 mas of the Herbig Ae star, 
and all but planetary mass companions exterior to O'.'S. We reach 5-0" contrasts 
limiting companions to planetary masses, 3-4 Mj at I'.'O and 2 Mj at 1'.'55 using 
the COND models. Collectively, these data strengthen the case for MWC 758 
already being a young planetary system. 

Subject headings: circumstellar matter instrumentation: high angular resolution 
polarization planetary systems: protoplanetary disks stars: individual (MWC 758) 
waves 
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Introduction 



We live in a planetary system with four giant planets, and over the past 15 years 
a wealth of other massive planets have been identified from radial velocity, transit, and 
microlensing studies, as well as from direct imaging of nearby stars. These studies have 
told us a great deal about the frequency of giant planets around older. Main Sequence 
stars, but do not provide insight into when, where, and how frequently they form in their 
natal circumstellar disks, beyond the trivial constraint that they must have formed while 
the disks retained abundant gas. A census of signatures of giant planet presence in young 
disks can be used to bound both the time required to form such bodies, and the portions 
of the disk occupied by giant planets. Such data are needed to constrain mechanisms for 
their formation and migration. Interferometric detection of young, Jovian -mass exoplanets 



and 



3rown dwarfs wit 



2012 



Huelamo. et al. 



lin t 



leir disks is challenging, but an emerging field (iKraus fc Ireland 
20111 ). however only a small portion of the disk is sampled. Search 
techniques focussing on the macroscopic, albeit indirect, signatures that such planets induce 
in their host disks offer an alternate way of carrying out such a giant planet census. 

Recent, high-contrast polarimetric differential imagery of the Herbig F star SAO 206462 
has demonstrated one such indirect si gnature : the prese nce of spiral arms interpreted as 



being excited by spiral density waves ( iMuto et al. 



20121 ). For SAO 206462, 2 arms were 



clearly resolved from the disk, and have amplitudes consistent with perturbe rs that are 
Saturn-mass objects. SAO 206462 is estimated to have an age of 9±2 Mvr (IMiiller et al. 



20091 ) and 



20111 ). older than age estimates for the formation of Saturn ( ICastillo-Rogez et al.l l 
potentially Jupiter in our Solar System. While tantalizing, these data, by themselves, do 
not provide insight into how common spiral features are in the disks of young stars, or at 
what point in the evolution of a disk they become detectable. 



The Herbig Ae star MWC 758 (HD 36112, A8Ve; V=8.27, B-V=0.3; H=6.56 



- 7- 



Beskrovnaya et al.lll999l ). (ASIVe lMeeus et al.ll2012l ) has also been identified as hosting a 



partially cleared cavity de tected in the sub-millimeter within 0'.'2 of t he star (llsel 



2010 



Andrews et al. 



201 ll ). The star has been dated to 3.5±2 Myr ( iMeeus et al 



a et al. 



2OI2I) 



comparable to the era of gi ant planet format i on in our So 



molecular and atomic gas fllsella et al 



201 



Salvk et al. 



accretion continues onto the star fiBeskrovnava et al. 



data place the star at 279 



nc fl 

-58 y. 



van Leeuwen 



the older Hipparcos measurement of 200^40 P*^ ( Ivan den Ancker et al 



ar System. The disk sti 



2011 



Meeus et al 



1 contains 



2OI2I), while 



9991 ). The revised Hipparcos parallax 



2007D. but much of the lite rature uses 



19981 ). To ensure 



compatibility with the older literature, we use both distance estimat es in this paper. Like 



SAO 206462, the disk has a low inclination from pole-on (i=21±2° 



Isella et al. 



2OI0I), 



facilitating detection of structure in the inner parts of the system. MWC 7 58 has a relative ly 



smal 



depletion of mm-sized dust grains in the inner portions of its disk (lAndrews et al. 



20111 ). and also has not been previousl y reported to host the distinctive mid-IR dip in the 



IR spectral energy distributi on (SEP) (llsella et a 



gas up to 2'.'6 from the star (jChapillon et al. 



20101 ). Millimeter interferometry traces 



20081 ). while millimeter dust emission can be 



traced ~1" from the star. The scale of the disk, the gap size, and the inclination of the 
outer disk all make this a suitable system to test the hypothesis that SAO 206462 is not 
unique in hosting spiral arms potentially associated with Jovian-mass planets. 

We report the first successful imaging of the disk of MWC 758 in scattered light, the 
discovery of spiral arms, the detection of a mid-IR dip in the IR SED as well as apparent 
mid-IR photometric variability, and constrain the mass of companions in and immediately 
exterior to the disk using a mixture of NIR high-contrast imagery, sparse-aperture masking 
interferometry at L', and the measured properties of the spiral arms. 
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2. Observations and Data Reduction 

In this section we discuss observations and data reduction for high-contrast imaging 
carried out using two instruments at the Subaru Telescope, a re-analysis of HST 
coronagraphic imagery, aperture polarimetry and photometry, and assembly of IR 
spectral energy distribution data, and VLT sparse aperture masking interferometry. The 
high-contrast Subaru imagery is processed both for optimal disk detection and for giant 
planet searches. 



2.1. Subaru/HiCIAO Polarimetric Differential Imagery 



MWC 758 was observed in the H-band 



imag ing instrument HiCIAO ( iTamura et al 



(X=1 . 635ta]f^ um) using t 



2006 



Hodapp et al 



2008 



le high-contra st 



Suzuki, et al 



20101 ) on the Subaru Telescope on 2011 Dec. 24 UT as p art of the Stra tegic Exploration 



of Exoplanets and Disk s with Subaru 



adaptive optics system ( iMinowa et al 



SEE DS) program (iTamura 



20091). The AO 188 



20101) provided a stable stellar point spread function 



fPSF, FWHM 



( iMarois et al. 



0'.'07). We observed with a combined angular differential imaging (ADI) 



20061 ) and polarization differential imaging (PDI) mode with a field of view 



of 10" by 20" and pixel scale of 9.5 mas pixel"^. A 0"3 -diameter circular occulting mask 
was used to suppress the bright stellar halo. Half-wave plates were placed at 4 angular 
positions from 0°, 45°, 22.5° and 67.5° in sequence with one 30 sec exposure per wave plate 
position for a total of 15 datasets, and a total integration time for the polarized intensity 
(PI) image of 1,440 s, after removing 3 low-quality images with FWHM> O'.'l, by careful 
inspections of the stellar PSF. The field rotation in this dataset was 9° on the sky, too 
small for ADI reduction of the data. While a PSF star was observed immediately following 
the PDI+ADI image set, the PSF had changed significantly from that observed during the 
observation of MWC 758, precluding conventional PSF subtraction. H-band photometry 
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for MWC 758, taken with the telescope immediately before the PDI+ADI dataset yielded 
H = 6.™'48±0.™01 in t he MKQ filter sys t em. T he polarimetric data were reduced using the 
procedure described in 



Hashimoto et al 



(120111) using IRA: 



2.2. Subaru/HiCIAO ADI observations 

A second HiCIAO dataset was obtained as part of the SEEDS program at 
(A=2.150±0.160 yum) on 2011 Dec. 26 UT, with 17 minutes of integration time in ADI 
mode. The field of view for this dataset was 20" x 20". The data were obtained in direct 
imaging mode, with no coronagraph, but with short 1.5s frame times to minimize image 
saturation. Only the peak of the PSF of MWC 758 was saturated in the data. The 
1.5s integrations were coadded into frames of typically 20 integrations, but after culling, 
some 10-integration frames were also included. After culling of bad frames, 640s of 
integration time were available. The data were processed using variants of Locally Optimized 
Combination of Images (LOCI) in two ways: using aggressive LOCI (na = 300, nfwhm = 
0.5) for a companion search and conservative LOCI (na =10,000, nfwhm = 1.5) to search 
for structure in the disk. A discussion of these approaches is in Thalmann et al. (2010). 
Either LOCI variant suppresses azimuthally-symmetric or large angle structure, resulting 
in voids compared to either the PI data or PSF-subtracted imagery. 

^ IRAF is distributed by the National Optical Astronomy Observatories, which are oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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2.3. Subaru/IRCS K' Data 



MWC 758 and a PSF reference star (HD 242067) were imaged on November 10, 
2011 with the Subaru Telescope using the Infrared Camera and Spectrograph (IRCS; 



Tokunaga et aL 



19981 ) and A0188 in natural guide star mode as part of program SllB-012. 
The data were obtained in the Mauna Kea K' filter (A = 2.12 fiin). The native pixel scale 
is 20.57 mas/pixel; we used the 0.15" diameter coronagraphic mask to reduce light from the 
PSF core spilling out to angular separations of interest for disk detection. Despite variable 
weather conditions, AO-188 delivered corrected images with a full-width half-maximum of 
0'.'07. 



Our data consist of coadded 15-second exposures taken in ADI mode (IMarois et al. 



20061 ) through transit HA = [-0.6,0.9] with a total field rotation of 127°. Our cumulative 



integration time for MWC 758 and HD 242067 is 38 minutes and 10 minutes, respectively. 



Currie et al. 



fcoilaU bh for IRCS, including 



Basic image processing follows steps outlined in 
dark subtraction, flatfielding, and distortion correction (corrected pixel scale = 20.53 
mas/pixel). We considered all HD 202067 frames to construct a reference PSF. 

Disk Imaging: To extract a detection of the MWC 758 disk, we explored two 
approaches: one using PSF subtraction of MWC 758 frames by a median HD 242067 frame 



and one leveraging upon an "adaptive" LOCI approach (A-LOCI, iCurrie et al. 



20m for 



this subtraction. In both cases, we first identified the highest-quality images of MWC 758. 
Specifically, we monitored the brightness of a background star 2" from MWC 758 (~ 15 
magnitude) and selected only the frames where the star's brightness is within 20% of the 
peak brightness of any frame. This selection criterion identified 50 MWC 758 frames (750 s 
total) to focus on for PSF subtraction. Further reduction steps are as follows: 



• Classical PSF Subtraction - Here, we subtract the median-combined HD 242067 image 
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from each MWC 758 image. We choose a weighting that minimizes the residuals 
between separations of 50 pixels and 150 pixels (~ l"-3"), though weightings over 
slightly different ranges in angular separations (i.e. 0'.'75-2'.'5) achieve similar results. 
We also performed a second subtraction weighting the reference PSF by an additional 
10% above the value that minimized the residuals for each MWC 758 image. We then 
derotated and median-combined each PSF-subtracted MWC 758 image. 

A-LOCI - Here, we treat all HD 242067 images separately as a "reference PSF 
library" . Instead of subtracting each reference image from the MWC 758 image at 
once, we follow the (A)-LOCI app roach, subdividing each science image and reference 



image into annular sections (see 



Currie et al. 



2012 



Lafreniere et al. 



20071 ) and 



determining coefficients for each reference PSF section that minimize the subtraction 
residuals within that section for a given science image. To better ensure that any 
"disk" signal is not instead due to PSF mismatch, we compute the cross-correlation 
function, Tcorr, for each pair of science-reference images and remove poorly correlated 
image pairs (we find the best results for rcorr > 0.85). We set the optimization area 
(in units of the image FWHM) over which we define the coefficients to A^a = 3000 
and adopt a radial width for the subtraction annulus of dr = 140, equal roug hly to 



the angular extent of the disk as seen in the submm (e.g. 



Andrews et all 120 11 



A-LOCI Limits on Planetary Companions - To search for companions to MWC 
758, we employed much more aggressive PSF subtraction settings using A-LOCI on the full 
set of images (150 images; 37.5 mi nutes). We adopted A-LOCI settings of 6 = 0.70, Na = 



250, dr = 5, and 



0.95 (see 



Currie et al. 



2012 



Lafreniere et al. 



20071 ). Additionally, 



^We detect the inner disk spirals as long as we set dr to be larger than the submm-resolved 
disk gap {dr > 25) 
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we employed a moving-pixel mask over the subtraction zone, using only the pixels outside 
of this zone to determine the LOCI coefficients to reconstruct a reference PSF annular 
region. We then quantify and correct for self-subtraction inherent in A-LOCI by comparing 
the flux of fake point sources prior to and after processing. We use the measured brightness 
of the star visible at 2" {K' ~ 16.93) and assume a brightness for MWC 758 of m.{K') = 
5.804 for flux calibration and thus to determine the planet-to-star contrast. 



2.4. Archival HST Coronagraphic Imagery 



MWC 758 was observed in the broadband optical twice by HST/STIS on 2000 Jan. 
16 and Jan. 19 as part of HST-GO-8474, with the star placed at a point along the STIS 
corongraphic wedge structure (wedgeAl.O) where the wedge occulted the inner 0'.'5. The 



STIS data (observation ids=O5KQ3010-20 and 05KQO 



the available library of PSF template data (iGrady et al 



10-20 ) were a poor color match to 



20051). vielding disk non -detections 



20051), a null 



with large color mis-match errors. When roll-differenced ( iLowrance et al. 
detection was also seen, indicating either nebulosity which was azimuthally-symmetric over 
15°, or concentrated within 0'.'5 of the star. The STIS data indicate only 3 point sources in 
the 25"x 50" field other than MWC 758, with the nearest at a separation of 2'.'1 at PA=309°. 

MWC 758 was also observed by HST/NICMOS as part of HST-GO-10177 on 2005 Jan. 
7. The observations consisted of two independent target acquisitions, short direct images 
in F171M, the coronagraphic observations at FllOW (Acerat=l-l=t0.3 /xm) with the region 
at r<0'.'3 occulted, and short direct-light FllOW imagery with the star unocculted before 
rolling HST by 29.9° and repeating the sequence. The direct-light imagery yielded FllOW 
magnitudes of 7.™314±0.014. The observations obtained at a given spacecraft orientation 
constitute a "vi s it". T he coronagraphic data for each visit were reduced as described in 



Schneider et al 



( 120061 ). and then a suite of PSF template data from the same HST cycle 
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were scaled, registered, and subtracted from the MWC 758 imagery. For the first visit 
for MWC 758 (visit 31), PSF template data from visits 3D (HD 15745, visit 2), visit 62 
and 61 (HD 83870), visit 53 (HD 35841 visit 1), visit 41 (HD 22128, visit 1), and 6B (HD 
204366, visit 1) provided good color and HST residual matches to MWC 758. For the 
second visit (visit 32), data from visit 48 (HD 38207 visit 1), 4C (HD 30447, visit 1),4D 
and 4E (HD 72390, visit land 2), 42 (HD 22128 visit 2), 62 (HD 83870 visit 2), 65 and 66 
(HD 164249, visits 1 and 2) , and 6B (HD 204366, visit 1) provided the best matches. The 
use of independent target acquisitions and largely independent suites of PSF template data 
mean that medians of the net images for each visit provide two essentially independent 
observation sets for MWC 758. Both final images provide no indication of nebulosity for 
r>0'.'8, with the bulk of the detected signal within 0'/5 of MWC 758. 

The angular scale of the nebulosity is too small compared to the NICMOS FllOW 
resolution element to assess azimuthal symmetry, but we have carried out aperture 
photometry between 0'.'3<r<0'.'5 for each of the final net visit images. The total brightness 
in those annuli in instrumental counts/s/pixel are 5183 for visit 31 and 5419 for visit 32, 
yielding only a ±2% peak to peak difference about a mean of 5301 counts/s/pixel. Using 
the NICMOS photometric calibration of lc/s/pixel=1.26xl0-^ Jy, the total FllOW flux 
density in the measurement annulus is 6.7 mJy. The 2MASS J magnitude for MWC 758 is 
7.22 (2.07Jy), so the total 1.1/xm light scattered by this portion of the disk is 0.32%. The 
nearby point source noted in the 2000 STIS imagery is recovered at r=2'.'215, PA=310.86°. 
It is not co-moving with MWC 758, and thus is a background object. 
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2.5. Infrared Survey Facility /Simultaneous IR polarimeter (SIRPOL) 

Aperture Polarimetry and Photometry 



NIR aperture p hotometry and 



2012 using SIRPOL flKandori et al. 



jolar imetry was carried out for MWC 758 on 8 March 



|2006[ ) at the IRSF. The SIRPOL observations were not 
accompanied by observation of a polarimetric standard star, so the 2MASS system was used 
to calibrate the NIR magnitudes of MWC 758 using field stars. We measured J=7.™1.93, 



H=6."560, Ks=5.'"804 with polarizations of 1.63% at J, 1.60% at H and 0.81% at Ks. T 



ne 



polarization levels are typical of low-inclination young stellar objects ( iPereyra et al. 



20091). 



2.6. The IR Spectral Energy Distribution of MWC 758 



MWC 758 is a well-studied Herbig Ae star with data from the far-UV through the 
millimeter. We have assembled a composite SED including lUE spectra SWP 53939 and 



LWP 30023, optical through M-band photomet ry flMalfait et al 



1998 



de Winter et al 



2001 



Beskrovnaya et al 



Bogaert et al. 



1994 : 



19991 ) ■ 2MASS, MSX, IRAS, an d ISO SWS dat a, 



2 epo chs of Broadband Array Spectrograph Sy stem spectrophotorn etry (IHackwell et al. 



Juhasz et al. 



mm , AKARI, and WISE 



19901), Spitzer IRS data originally published by 
data (excepting band 2, 4.5 /xm). We also include the HST/NICMOS FllOW photometry 
as well as the photometry obtained in tandem with the HiCIAO PDI observation and the 
SIRPOL data. The data, color-coded by source, are shown in Figure 1, together with 
a Kurucz T=8250 K photospheric model, reddened by an R=3.1 extinction curve with 



E(B-V)=0.13. A similarly good fit at short w avelengths can 



reddening and with Te//=7580K, as used by 



3 e ach ieved without foreground 



Andrews et al. 



( I2OIII ). This model requires a 



contributio n from UV excess emissi on (e.g. accretion luminosity) to match the observed 



FUV fiux (IMartin-Zaidi et al. 



20081), but is consistent with low line-of-sight extinction 



(E(B-V)<0.1) for the source to be detected by FUSE. 
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2.7. VLT/NACO Sparse Aperture Masking Interferometry 



MWC 758 was observed on 11 March 2012 with VLT/NACO Sparse- Aperture 
Masking^. Observations were taken in the L' band (Al'= 3.80 ± 0.31 fim) using the "7 
holes" aperture mask and the IR wavefront sensor (WFS). A total of 24 data cubes, each 
one made of 330 frames of 0.12s integration time were obtained for MWC 758 in L' on the 
night of 2012-03-12(UT). Observations of two cahbrator stars , HR 1921 and HD 246369, 
were obtained for a total of 8 data cubes, interleaved every 8 MWC 758 datacubes. The 
MWC 758 observatio ns were processed with the Observatoire de Paris SAMP pipelines 



(ILacour et al 



2011a 



The use of the "7 holes" (C7-892, 



Tuthill et al. 



20101 ) aperture mask transforms the 



telescope into a Fizeau interferometer. The point spread function is a complex superposition 
of fringes at given spatial frequencies. In specific cases, p upil-masking can outperform mor e 



traditional differential imaging for a number of reasons (ITuthill et al. 



2006 



Lacour et al. 



2011al ). First, the masks are designed to have nonredundant array configurations that 
permit phase deconvolution; slowly moving optical aberrations not corrected by the AO 
can be accurately calibrated. Second, the mask primarily rejects baselines with low spatial 
frequency and passes proportionately far more baselines with higher A/B (where B is 
baseline length) resolution than does an orthodox fully filled pupil. Third, high-fidelity 
recovery of phase information allows "super resolution" , with a marginal loss of dynamic 
range up to A/2D (where D is the mirror diameter). The principal drawback is a loss in 
throughput so that photon and detector noise can affect the signal-to-noise ratio even where 
targets are reasonably bright for the AO system. The effective field-of-view of SAM is 
determined by the shortest baseline so that the technique is not competitive at separations 
that are greater than several times the formal diffraction limit. For more details on the 
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SAM mode, please see e.g. 



Lacour et al 



(1201 Ibh 



Tuthill et al 



mm . 



Results 



In this section we discuss the disk detection, compare it with synthetic imagery for 
a model fitting the sub- millimeter conti nuum data and the IR SED, report the discovery 
of spiral arms in the disk, and following 



Muto et al. 



( I2OI2I ) discuss the implications for a 



giant planet capable of exciting the more completely imaged spiral arm. We next discuss 
limits on stellar, brown dwarf, and giant planetary mass companions provided by the sparse 
aperture masking data and the Subaru observations. 



3.1. MWC 758 as a Gapped Disk 



Previous studies of MWC 758 had identified this system as hosting a cavity visible 
in sub-millimeter continuum interferometry, but lacking a very conspicuous dip in the IR 



SED near lOyum ( llsella et al. 



201 



Andrews et al. 



2OIII ). In contrast, our composite SED 



based on an assembly of archival photometry and spectrophotometry clearly shows a dip 
(fig. 1), with indications of mid-IR variability, providing the first independent confirmation 
that MWC 758 is what has been termed a gapped or "pre-transitional" disk with an inner 



disk structure (Isella et al. 



(llsella et al. 



20081), a wide 



gap, and an outer d isk seen in millimeter data 



2OIOI ). and thermal emission (IMariiias et al.ll201ll ). When compared with other 



transitional disks, the depth of the 10 fim dip is smaller, consistent with the lack of a cavity 
in our FllOW, H, K' and imagery. Modeling of other disks with HiCIAO observations 
have suggested that the absence of a cavi ty at H is not at ypical of transitional disks, and 



can be accounted for by grain filtration (IRice et al. 



20061 ). with large mm-sized grains 



deeply depleted in the cavity region, while smaller grains, which are more tightly coupled 
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to the K less depleted flDong et al.ll2012l : IZhu et al.ll2012[ ) 



3.2. The Structured Disk of MWC 758 

The disk of MWC 758 is detected in polarized, scattered light from 0'.'2 to O'.'S (223 AU 
assuming d=279 pc or 160 AU assuming d=200 pc, fig. 2) in polarized inte nsity at H-band 



(fig. 2). The disk detection did not require unusually favorable conditions (IKusakabe et al. 



2012 



, e.g. MWC 480, ): H-band photometry was in the range of previous NIR data for the 
star (fig. 1). The detected nebulosity is highly asymmetric about the star, with the greatest 
extent seen to the west of the star at H, while to the east, nebulosity is not detected exterior 
to O'.'S. To the west of the star, the scattere d light disk extends as far as the sub-millimeter 



continuum can be traced (jlsella et al 



20101) (fig. 3a,b). The disk is also detected at (fig. 



2c) to within O'.'l of the star, and in K' (fig. 2d) to within 0'.'25, and at l.lyum (fig. 2e). A 
false-color composite of the H PI (blue) and K' (red) is shown in figure 2f. This composite 
image also allows us to identify a region of reduced polarization intensity along PA=155±5° 
and PA=330±5°, aligned with the projection of the disk semi-minor axis. To the east of 
the star, the H PI surface brightness ocR"^'^^"'^ (fig. 3). To the west of the star, it drops as 
Qj.j^-2.8±o.i^ consistent with the disk exhibiting some fiaring. 



Unlike other protoplanetary ( IKusakabe et al.ll2012l ) or transitional disks ( iMuto et al 



2OI2I ) imaged with HiCIAO in H PI, the refiection nebulosity is not aligned with the disk 



semi-major axis (PA=65±7°). Moreover, we detect scatte red light within t he region of 



20101). rather than 



partial clearing reported in the sub-millimeter continuum ( jlsella et al.l 
the expected cavity. The PI for r>0'.'2 is 6.095±0.085 mJy, or 0.09% relative to the star. 
Between 0'.'3 and O'.'S the PI is 2.72 mJy, or 0.0425% relative to the star, indicating that 
the region of the arms is providing half of the observed polarized intensity. If we assume 
neutral scattering, the PI in the arm region 0'.'3<r<0'.'5 can be compared with the FllOW 
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total intensity of 0.32%, indicating a polarized light fraction of ~13.3% which is in the 
range of polarization fractions observed for other Herbig Ae stars. 



3.3. Comparison with Monte-Carlo Radiative Transfer Models for the Disk 



Given the structure seen in the H-band PI imagery, it is natural to explore what 
we would have expected to have seen based on combined modeling of the sub-millimeter 
data and the IR SED. To address this goal, we carried out M o nte Ca rlo radiative transfer 



modeling using the code develope d by 



as described in 



Dong et al. 



Whitney et al. 



fl2003a 



2OI2I ). with modifications 



(120121 ). In particular, we adopt a disk ou ter radius of 200 AU 



viewed at i=23°, and assume an accretion rate of M = 10 ^Mq yr ^ (lAndrews et al. 



mm. 



We model the entire disk with two components: a thick disk with small grains (~ /xm-sized 
and smaller), which represents the pristine dust grains from the star forming environment, 
and a thin disk with large grains (up to ~mm-sized), which is the result of dust growth 
and settling towa rd the disk mid plane. We adopt the interstellar m edium dust model in 



Kim et al. 



(119941 ) for our small-grain dust, and the dust model 2 from 



Wood et al. 



mm for 



our large-grain dust, which has a power law size distribution up to 1 mm (with 
index 3.5). The large-to-small-dust mass ratio is ass umed to be . 85/0. 15 as in 



power law 



(120121 ). The properties of the grains can be found in 



Wood et al 



(12002). 



Dong et al. 



With these assumptions, the total dust mass of the disk is 0.067Mj. The radial surface 
density distribution of both dust grain populations is assumed to have the following profile 



E(i?) 



/lAuy 



where Sq is the normalization factor, R is radius, and a is the power law index. In the 
vertical direction, we assume 



11 



(2) 
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where p{R, z) is the local volume density, z is the vertical dimension, and H is the scale 
height (input parameters in the code). Radially, both scale heights for the small and large 
dust vary with radius as 

H oc i?^, (3) 

where /3 is a constant power law index, assumed to be 1.08 in this work. The scale height, 
h, of the small dust is assumed to be 15 AU at 100 AU, while that for the large dust is 
assumed to be 1/5 of the h for the small dust, to reflect the fact that big grains tend to 



settle to the disk mid-plane, while sm all grains tend to be we . 



have a much larger vertical extension (IDullemond fc Dominik 



1 coup l ed with the gas and 



20061 ). 



2004b, 



2005 



D'Alessio et al. 



Our disk model has a gap structure located from 0.5 AU to 73 AU. Inside the gap 

ly resolved SMA obs ervations, 



20111 1. However, 



there are no large dust grains, as suggested by the spatial 
which revealed a clear cavity 73 AU in radius at 880 /im (j Andrews et al.l 
the H-band PI data demonstrate that there is a significant amount of scattered light within 
this radius, indicating the presence of starlight-scattering small dust grains. We set a = 1 
for both dust populations in the outer disk, and a = —3 for the small dust inside the gap, 
wi th a continuous su rface density distribution across the gap edge. As discussed in detail 
in 



Dong et al. 



( l2012l ). this kind of surface density distribution of the small dust results in 
a heavy depletion at the inner disk, producing the IR excess depression around 10 fim on 
the SED seen in transitional disks, but a smooth NIR scattered light image with no breaks 
or discontinuities of the surface brightness radial profile at the inner edge of the outer 



20081) 



disk. MWC 758 has a NIR excess (fig. 1), and VLTI/AMBER data fllsella et al. 
indicating a small amount of hot, small dust grains at the inner-most radii. To reproduce 
these features, we put a inner rim located at 0.33 to 0.5 AU. The surface density of the rim 
is ~ 0.002g/cm2. 



-20- 



Figure 4 shows the model SED and iJ-band polarized light image from the modeling, 
with the disk major axis horizontal, and the near side of the disk in the bottom half of 
the image. The model ima ges show both the PI as modeled, and after convolution with an 



observed HiCIAO PSF(see 



Dong et al. 



20121 for details). The model reproduces the extent 
of the polarized intensity to the west of the star, but, as expected, predicts a similar extent 
of the scattered-light disk to the east of the star, beyond O'.'S. The model predicts a bright 
arc on the near-side of the image, which is not seen in the data. The raw model imagery 
predicts a cavity visible at H in the immediate vicinity of the star, which is not conspicuous 
after convolution with the PSF, and which is in the region occulted by the coronagraph in 
our H-band data. Higher Strehl ratio and smaller inner working angle imagery, such as can 



be provided by extreme AO systems, will be required to test these predictions of the model. 



3.4. Spiral Arms 



The most distinctive features in the NIR images which are not captured in the model 
imagery are a spiral feature on the east side of the star, wrapping toward the S, which we 
term the SE arm, and a similar feature originating in the NW and wrapping to the N, 
which we term the NW arm. The spiral arms are most clearly detected interior to O'.'S (100 
AU for d=200 pc, or 140 AU for d= 279 pc). Our data for MWC 758 therefore demonstrate 



that SAO 206462 flMuto et al 



2OI2I ) is not unique in hosting spiral arms straddling the 



region o 



the partially c 



142527 (Casassus et al. 



eared millimeter cavity. M oreover, like SAO 206462, and unlike HD 



2012 



Rameau et al. 



20121 ) there is no indication of a cleared zone 



at H-band for r>0'.'2 (40 AU for d=200 pc, 56 AU for d=279 AU). 

The spiral arms in the disk of MWC 758 are detected as a contrast in surface 
brightness. The drop in polarized surface brightness immediately exterior to the SE arm 
is consistent with that arm fully shadowing the outer disk. Shadowing by the NW arm is 
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less complete: some signal is seen beyond the arm in our 2011 data. We measure the SB 
contrast of the arm features by measuring the amplitude at the SE arm at PA=90° and the 
corresponding surface brightness the same distance from the star at PA=270°. For the SE 
arm we find a ratio of 50%±20%. The angular extent of the NW arm is sufficiently small 
that to confidently measure the amplitude and fit the pattern we will require data obtained 
with a combination of higher Strehl ratio, and a coronagraph with smaller inner working 
angle, such as may be provide d by SCExAO in tandem with A0188+HiCIAO on Subaru 



(IMartinache et al 



2011 



2012I ). We therefore restrict all but qualitative discussion of the 



arms to the SE arm. 



Both of the arms have the same rotation sense. If they are trailing structures, as 
expected for spiral density waves, the rotation of the disk, and any companions asso ciated 



2010h . then 



with the arms, is clockwise. Comparison with the millimeter data (llsella et al.l 
indicates that the northwest side of the disk is the near side of the disk. Our PI image 
is fully consistent with the disk being viewed at a low inclination from face-on, since it 
lacks t he compression of the PI toward the disk semi-major axis which is conspicuous at 



i>38° ( Kusakabe et al. 



2OI2I ). Given the asymmetric i llumination of the disk, we adopt the 



millimeter disk inclination i=21±2° (jlsella et al. 



2010h . 



3.5. The disk at other wavelengths 

To be correctly interpreted as spiral density waves, the spiral arms should be recovered 
in observations made at different wavelengths and with different instruments. The IRCS K' 
data detect the disk in scattered light to at least O'.'S from the star, and recover both spiral 
arms using either classical PSF subtraction and using A-LOCI as long as we set dr to be 
larger than the size of the disk gap {dr > 25). The data sample the disk less completely, 
due to the very short integration times for each individual exposure. However, they recover 
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the inner portions of the arms, in the conservative LOCI processing of the data, to an inner 
working angle of O'.'l but suppress them in aggressive LOCI processing. 

The brightest portion of the disk as seen at 1.1 /im is the region occupied by the spiral 
arms, between C'S and O'.'S. The HST/NICMOS data from 2005 provide no indication of 
the asymmetric illumination seen in the 2011 data: this may either be due to changes in 
illumination over time, or to contrast limits of the HST/NICMOS data. However, we note 
that asymmetric illumination of the disk was not seen in the HST/STIS data from 2000. 
Additional observations will be required to establish whether the illumination of the outer 



disk changes with time. The locations o: 



from the MWC 758 disk flMarinas et al. 



the arms coincide with the brightest 12/im signal 



20111 ). The coincidence of structural features in 



scattered light, polarized scattered light, and in thermal emission indicate that the spiral 
arms are density features, and not merely distortions of the disk surface. The rapid drop in 
surface brightness seen in the H-band PI data to the east of the spiral arms (fig. 3) further 
indicates that the arm is optically thick and extends sufficiently far above the undistorted 
disk surface to cast a long shadow. The reg ion of the NW arrn near PA=340° coincides with 



the clump in continuum emission noted by 



Isella et al. 



(l2010l ). although with a synthesized 



beam size of 0'.'76x0'.'56, the available sub-millimeter continuum interferometry lacks the 
angular resolution to directly detect the arms. As the angular resolution of sub-millimeter 
interferometry improves, we expect that these features should be detectable in continuum 
and gas observations, such as can be provided once ALMA reaches its full extent. 



3.6. Fitting the Spiral Arms 



We follow the approach given in iMuto et al.l (120121 ) to fit the spiral arms. Our initial 
fits have been restricted to the H-band PI data. Given the nearly face-on system inclination, 
we have not de-projected the data to compensate for the non-zero inclination, but have 
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compensated for the drop in illumination of the disk (fig. 5). From the data, we have 
looked for the local maxima of the surface brightness xr^ profile for each radial profile 
both for the SE spiral and the NW spiral. If the features are a part of a ring, the points 
would lie horizontally in a polar coordinate plot of the intensity. The SE spiral is clearly 
non-axisymmetric, but only a small angular extent of the NW arm is non-axisymmetric, 
given the inner working angle of our data. For the remainder of this study we restrict our 
attention to the SE spiral. 

Assuming i=21° and PA=65° for the disk major axis, the fitting function for the spiral 

is: 



e{r) 



sgn(r — Tc 
he 



X 



-7 



(4) 



1 + 5 1-7 + 5 V^c/ j \l + 5 1-7 + 5 

Here, as in SAO 206462 paper. The corotation point (~ launching point of the spiral) 
is at {rc,6o)- The disk rotation profile is Q{r) oc r~"' and the disk sound speed profile is 
c(r) oc r~^. The disk aspect ratio (H/R) at the corotation radius is he- Note that assuming 
H = c/Q, (3 that appears in Equation (3) is related to 7 and 5 by /3 = 7 — 5. We also 
note that if the disk temperature varies as T oc r~'^, q and S are related hj 6 = q/2 since 
temperature is proportional to the square of the sound speed and the disk sound speed 
profile is c(r) oc 

Note th at the sign a f ter th e first ^0 in the right hand side of equation (4) is different 



( I2OI2I ). due to the counter clockwise rotation of the SAO 206462 



from that in iMuto et al. 
disk. 

There are five fitting parameters: {rc,6o,hc,'y,S). Among them, three {hc,'~f,S) are 
determined by the disk structure and the two (rc, 60) determine the location of the spiral 
feature. Assuming Keplerian rotation, 7 = 1.5. In our study of SAO 206462 (Muto et al. 
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2012), the sound speed profile S was fixed, but we have varied this parameter as well. There 
is a range of values related to 6 in t he literature, ranging; between 0.45 (the parameter 



in Andrews et al. 2 011), 0.19 (q in 



m 



Isella et al. 



Chapillon et al. 



20081 ) . and 0.05 (the ( parameter 



20101 ). We note that for 5=0.5, the opening angle of the disk is almost 
constant, in disagreement with our radial SB profile. We looked for the best-fit parameters 
in the least sense by dividing the (rc, ^o) space in 128 x 128 cells and {he, S) space in 
100 X 100 cells with equal spacing in the linear scale. The parameter space explored, and 
best-fit parameters, a re shown in table 1. Our value for h^ is similar to that derived by 



Andrews et al. 



( I2OIII ). The external perturber fit resulted in a reduced of 0.68. 



Now we consider the degeneracy of the parameter space. The degeneracy of parameters 
of the spiral position (rc, ^0) and those of the disk parameter [he, 6) are investigated 
separately. The left panel of Figure 6 shows the degeneracy of (rc,6'o). The best fit to the 
SE spiral arm places the launching point of the spiral outside of both where we detect 
the spiral arm, and at a radius exterior to the dust disk as seen in the sub-millimeter. It 
might worth pointing out that there is a range of parameters where the corotation radius 
Tc is inside the spiral arms and where the values of becomes small, although these 
parameters are outside the range of 99% confidence. 

In the right panel of Figure 6, we show the degeneracy in the {he, S) parameter space. 
For each fixed values of {he, 5), we looked for the set of {tcOq), which gives the minimum 
5x^ value. Note that 5 = is a constant temperature, fiared disk and 5 = 0.5 is the 
constant opening angle disk. Although the degeneracy is significant, we favor a disk with a 
rather fiat, constant opening angle geometry. In this case, the spiral can cast a long shadow 
over the outer disk rather easily, if the spiral has sufficient amplitude, as is observed in the 
H-band PI imagery. We also note that the parameters with warmer disk {he > 0.1) are 
favored. 
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3.7. Estimating the Perturber Mass for the SE Spiral 

If we assume that a perturbing body excites the spiral arm, the mass of the perturber 
launching the spiral density waves can be estimated independent of knowledge of where it 
might be located, if there are data indicating the relative disk scale height {he), the mass of 
the star, and the amplitude of the spiral pattern as follows: 

Mp/M*= (Pattern Amplitude) {h^f 

Adopting 1.8±0.2 Mq for the stellar mass, /ic=0.18 from the pattern fitting for the SE 
spiral, and a pattern amplitude of 0.5±0.2 from the surface brightness contrast for the SE 
arm, we find a perturber mass for the SE pattern of ~5!!^4 Mj, consistent with the perturber 
being a giant planet, but excluding stellar or brown dwarf mass objects. This mass estimate 
is also consistent with accretion continuing onto the star. 



3.8. Contrast limits on Companions to MWC 758 within or near the disk 

The 5-0" contrast map from our sparse-aperture masking (SAM) NACO observations 
is presented in figure 7a. At L' band, SAM is sensitive to companions at separations <300 
mas from the primary. Attained contrasts range from Amag(L) = -4 to Amag(L) = -8. 
Contrasts we r e converted to minimum detectable companion mass using the models of 



Baraffe et al 



(11998 



2OO2I ). and adopting an age of 3.7 Myr and a distance of 279 pc. Similar 



limits are achieved for 5 Myr and a distance of 200 pc. With these assumptions, minimum 
detectable mass images are shown in figure 7b. Our SAM observations are sensitive down 
to very low mass stellar companions (~80 Mj): we conclude that MWC 758 is a single star. 
Thus, our limits for companions within the region occulted in the HiClAO H-band PI data 
are at the low end of the stellar/ high end of the brown dwarf range. 



Figure 8a shows S/N maps for the and K' datasets. Within 3" of the star, the 



- 26 - 



only point source object is the K=16.98 magnitude background star, which was detected 
at 37(7 at K^. K-band 5-cr contrast hmits for the data (from 0'.'l-0'/25) and the K'data 
(r>0'.'25) are shown in figure 9. The deepest imagery in our set of observations is the IRCS 
K' data which provide 5a contrasts of 1.4x10"^ at 0'.'25, 2.1x10'^ at O'.'S, 3.3x10"^ at 1", 
and 1.12x10-6 at 1'.'55 (see fig. 8). 

Adopting an age of 5 Myr and the COND models (Baraffe et al. 2003), these contrasts 
correspond to 5 a detection limits of 15-20 Mj at 0'.'25 , 8-10 Mj at O'.'S, and 3-4 Mj at 1", 
and near 2 Mj at 1'.'55 for distances of 200 and 279 pc respectively. These data exclude 
brown dwarf-mass companions exterior to O'.'S and are in general agreement with the mass 
estimates based on the spiral arm fitting. Interior to O'.'S, our upper limits are consistent 
with at most a low-mass brown dwarf close companions to MWC 7S8. Our mass limits are 
consistent with accretion conti nuing onto the star, as observed, and as pr edicted for disks 



with Jovian-mass companions ( iRice et al. 



2006 



Lubow fc D'Angelo 



20061 ) 



4. Discussion 



Confirmation that MWC 758 is a gapped disk: Historically, transitional disks have 
been defined by the weakness of the IR SED near lO/xm, compared to the excess at longer 
wavelengths (and in some systems at shorter wavelengths), in addition to the detection of a 
wide gap or cleared cavity in sub-millimeter continuum aperture synthesis imagery. Figure 



1 demonstrates the presence of a dip in t 



indicating an inner dust belt (ilsella et a 



in millimeter interferometry ( Ilsella et al. 



le IR SED near lO/xm. Together with VLTI data 



20081) . and the central cavit y and outer dis 



K seen 



20101 ) ■ and thermal emission ( iMarinas et al 



201ll ) 



the SED data indicate that MWC 7S8 is a transitional disk. When compared with other 
transitional disks, the depth of the 10 fim dip is smaller, consistent with less complete 
clearing of the few to 100 AU region. However, the depletion of millimeter-sized grains 
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flAndrews et al.ll2011f ) is similar to SAO 206462. 



Spiral arms are not unique to SAO 206462: MWC 758 is the third Herbig Ae star 
to show spiral arms i n the inner parts of the disk, extending from the region of the 



sub-mill imeter cavity 



142527 (ICasassus et al. 



sella et al. 



2012 



2010) into the outer disk. Like SAO 206462 and HD 



Rameau et al. 



20121 ) . the disk is dynamically warm (he >0.1) 



indicating that spiral arms should be detectable in other warm disks, if Jovian-mass planets 
are present in the disk, and the region of the (sub) millimeter cavity retains at least some 
small-grain material. We identify potentially 2 arms in the disk of MWC 758, although we 
fit only the one with the larger angular coverage. The spiral arms are seen H, K' and K^, 
and potentially also at 12//m, demonstrating that they are visible both in scattered light 
and thermal emission. While unresolved, they are also detectable in HST/NICMOS 1.1/im 
imagery. MWC 758 provides the first clear-cut case where an arm is optically thick, and 
shadows the outer disk, accounting for the disk non-detection with HST/STIS (IWA=0'.'5). 
The lack of similar shadowing of the west side of the disk suggests that the NW arm has a 
smaller amplitude, compatible with a lower mass. 

Limitations of sparse- aperture masking in companion searches for Herbig Ae stars: 
Sparse-aperture masking trades contrast for resolution, realistically obtaining contrasts of 
5-8 mag in the inner 0.1". This is sufficient to detect brown dwarf and planetary companions 



to yo ung, solar-like stars such as LkCa 15 and T Cha (iKraus fc Irelau' 



2012 



Huelamo. et al. 



201 ll ). However, at ages >3 Myr, contrasts of >10 mag are required to detect exoplanet 
and low mass (<30 Mj) brown dwarf companions to young, intermediate-mass stars such as 
MWC 758. Nonetheless, a close-in brown dwarf or ste llar companion can have a significant 



impact on disk properties (see, e.g.. 



Biller et al. 



20121 ). Thus, constraining the existence of 



such a companion is key for realistic disk modeling of the full system. 



Comparison of mass and luminosity constraints: Our mass estimate for the perturber 



- 28 - 



of the SE arm is consistent with the presence of at least one Jovian-mass, but sub-brown 
dwarf mass planet in this system. The width of the sub-millimeter cavity and the dip in t he 



IR SED are consistent with dynamical clearing by more than one body (IZhu et al 



2Qll|), 



as is the presence of an inner disk with an inclination which is 20° larger than for the outer 



disk (llsella et al. 



2008 



20101 1. The spiral arm appears to continue in to O'.'l, based on the 
Kg imagery, but characterizing the spiral requires imagery with higher Strehl ratio and 
smaller inner working angle than the data presented in this study. We estimate the mass of 
the SE arm perturber, based on surface brightness contrast, for the SE spiral to be 5^1 Mj, 
consistent with the pattern being launched by a Jovian-mass planet and excluding a stellar 
or brown dwarf-mass perturber, independent of whether the perturber is in the disk or 
exterior to it. This is broadly consistent with the SAM constraints excluding a stellar mass 
companion within O'.'S of the star, and ADI contrasts excluding brown dwarf mass bodies 
in the outer disk. Our 5cr contrast limit for an external companion is consistent with our 
perturber mass estimate, but does not inform on whether the perturber is has a luminosity 
compatible with the COND models, or is und er-luminous relative to th ose models, as 



Spiegel &: Burrows! (120121 ). Higher contrast 



predicted for cold-start models of giant planets 
imagery will be required to constrain the luminosity evolution of any planetary-mass bodies 
in the disk of MWC 758. 

Comparison with SAO 206462: MWC 758 resembles SAO 206462 in hosting 2 
spiral arms which are trailing with respect to the disk rotation, and which are broadly 
consistent with being associated with Jovian-mass perturbing bodies. If the spiral arms are 
independent patterns and are associated with giant planets, the presence of 2 arms each in 
SAO 206462 and MWC 758 and at least 4 in HP 142527 supports predictions that wide 



gaps must be cleared by multiple planets ( IZhu et al 



20111). Our fit to the SE spiral arm 



suggests a perturber mass ~10x larger than for the SAO 206462 perturbers, consistent with 
shadowing of part of the outer disk. If sufficiently many transitional disks with spiral arms 
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can be imaged, it will ultimately be possible to compare the giant planet mass function with 
studies of radial velocity giant planets, and to search for differences in the mass function as 
a function of stellar properties or system age. 

Even at this early stage in the search for giant planets in transitional disks, the 
prese nce of spiral arms in systems with ages ranging from 3.7-5 Myr (MWC 758) to 9±2 



Mvr (iMuUer et al. 



( lAndrews et al 



2011 



SAO206462) but otherwise similar levels of large grain depletion 
20111 ) places constraints on whether there is a characteristic age for giant 
planet formation and disk clearing. In particular, our data can be reconciled with theoretical 
predictions that once gaps a re opened, the inner disk should rapidly drain onto the star 



( lErcolano. fc Koepferl 



2OI2I ). if and only if giant planet formation is a process with either a 



variable onset and/or duration. 

A Planet and Inner Hole Formation: If the observed spiral structure is attributed to 
the existence of an unseen planet, such a planet might not account for the opening of the 
inner hole, since the location of the launching point is favored to be outside the location 
of the spirals (or the disk). Another possibility is that a planet which excites the spiral 
structures might actually be scattered from the inner disk, where the giant planet formation 
timescale is more rapid. In this case, the spiral fitting presented in this paper is relevant 
if the planet's eccentricity and inclination are well damped. If it still retains significant 



eccentricity, the shape of the spiral is mo re complicated and a. 



variability within one orbit of the planet f lCresswell et al. 



2007 



■ SO there is a strong t ime 



Bitsch fc Kiev 



2OIOI). In 



this case, the fitting model for the spiral should be modified. We also note, that a planet 
external to the disk, if confirmed and found to be on an orbit indicating scattering from the 
inner disk, would imply that the disk of MWC 758 must host additional giant planets. 

Implications for the future: MWC 758 has a steadily increasing body of data that 
indicates that the disk of this Herbig Ae star hosts at least 2 giant planets which are 
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dynamically sculpting the disk, both by clearing a wide partial gap in the disk, and by 
exciting spiral density waves which can be traced in our data to 100 -140 AU (depending on 
d=200 or 279 pc) from the star. Such waves are expected to pr ovide a means of funneling 



disk material onto and past any planets, until it reaches the star (IDodson-Robinson fc Salyk 



20111), consistent with both theoretic al predictions for the point at which a corn panion 



can choke off accretion onto the star fiRice et ali 12006 



presence of on-going accretion onto the star (iBeskrovnaya et al. 



Lubow fc D'Angelo 



19991). 



20061 ) and the 



Higher Strehl ratio/angular resolution imagery are required to map the NW spiral in 
this system. Deeper ADI observations are also needed to either directly detect planets or 
to place more stringent luminosity limits on them. Such advances should be possible as 
the next generation of AO systems and coronagraphs commission in the course of the next 
year. Further advances are expected in the sub-millimeter: a key assumption made in this 
study is that the surface brightness contrast of the spiral arms is a reliable proxy for density 
contrast measurements made at wavelengths where the dust disk is optically thin. Such 
data with sufficient resolution to separate the arms from the bulk disk are required, and 
can be obtained once ALMA reaches its full capabilities, in late 2013. 
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MWC 758 




X (/i-nn) 



Fig. 1. — Composite SED of MWC 758. The star has hmited UV data, optical data indicating 
minor photometric variabihty, consistent with a low inclination, NIR data showing a larger 
range, and mid-IR data from a variety of sources indicating a factor of 2-3 flux variation from 
5-30 Mm ■ The mid- IR data, particularly the Spitzer IRS spectrum previously presented by 



Juhasz et al 



(|2010[ ). clearly demonstrate the dip in the SED characteristic of pre-transitional. 



or gapped disks. 
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f) 



OS 




2 -0.8 -0.4 0.0 0.4 0.8 1.2 



b) 



Fig. 2.— The disk of MWC 758: a) the 880/im continuum after llsella et all (120101 1 . the 
dashed hne indicates the disk semi-major axis, b) H-band polarized intensity, c) data 
with conservative LOCI processing, d) K' intensity with conservative LOCI processing, e) 
HST/NICMOS total intensity data following PSF subtraction, f) color composite of H PI 
and K' data. The images are all 2'.'4 on a side, and oriented with north up and east to the 
left. The HiCIAO H PI data use an 0'.'3 diameter coronagraphic spot, the NICMOS data 
uses an 0'.'6 spot. The IRCS data use an O'.'IS diameter occulting spot, and the data have 
an inner working angle of 0.1". 
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Table 1. Spiral Arm Fitting 



Parameter Search Range Best Fit External Perturber 



rc 0'.'05 < rc < 1'.'55 = 1'.'55 

e 0<eo<2n 6*0 = 1.72[rad]^ 

he 0.05 <hc< 0.25 he = 0.182 

S -0.1 < (5 < 0.6 6 = 0.6 

^measured from PA=90° 

MWC 758 Rodial Profile 




-0.8 -0.6 -0.4 -0.2 0.0 0.2 

Log(r) orcsec 

Fig. 3. — Radial surface brightness profiles at PA=270° (black) and PA=90°(gray) together 
with power-law fits. 
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0.3 1.5 6.2 25.0 100.0 

Fig. 4. — Radiative transfer modeling results for MWC 758, showing the SED (left), the raw 
model if-band polarized light image (middle) with a cavity visible in the inner disk, and the 
convolved if-band polarized light image (right). The images are in unit of mJy/arcsec^. 
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Fig. 5. — Top left: the surface brightness of the scattered hght normahzed by r^, where 
r is the distance from the central star. The shape of the spiral is more enhanced, and this 
data is used for the analysis for the spiral fitting. Top right: the points that represent the 
spiral shape is shown by red points, which are used as data points for the fitting. Bottom 
left: the shape of the spiral given by Equation (4) with the best fit parameters (Table 1) is 
shown by green curves. The red points are the same with those in the top right panel. 
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Fig. 6. — Left: The degeneracy of (rc, ^o) parameters plotted in the 2D spherical coordinate. 
The color shows the values of 5x^ with respect to the best-fit values, which are shown by the 
red cross. The contour of 99% confidence level is shown by the green line, and the black circle 
is the radius of r = 1'.'55, which is the upper limit of the parameters that are searched. Right: 
The degeneracy of (/ic, 5) parameters. As in the left panel, the contour of 99% confidence 
level is shown by green curves and the best-fit parameter is indicated by the red cross. 
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Fig. 7. — left: 5-cr contrast map from our SAM NACO observations. Our observations are 
sensitive to companions 4 to 8 magnitudes fainter than the primary star at separation of 300 
mas or less, right: Minimum detectable companion mass m ap. Contrasts were conver ted to 



Baraffe et al. 



(11998 



2OO2I), and 



minimum detectable companion mass using the models of 
adopting an age of 3.7 Myr and a distance of 279 pc. We are sensitive to companions right 
down to the star / brown dwarf boundary, i.e. ~80 Mjup- From these observations, we 
conclude that MWC 758 does not have a low mass stellar companion within 300 mas of the 
primary star. 
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Fig. 8.— Aggressive LOCI S/N map for MWC 758 at K, (left) and K' (right). The field 
shown for the data is ~2'.'5, while that for the K' data has an outer radius of 3'.'04. The 
only point source detected in either dataset is the background object first seen by HST. 
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Fig. 9. — Five a contrast limits for Jovian- mass bodies in or near the disk of MWC 758 
derived from the K, data (0.1-0.25"), and the K' data (r>0.25"). The discontinuity at 0'.'25 
reflects the change in exposure depth between the and K' datasets. 



